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We theoretically examine the superconducting state of BaFe2(Asi- :E P a ;)2, an isovalent dop- 
ing 122-iron-pnictide superconductor. We construct a three-dimensional ten-orbital model by 
first-principles band calculation, and investigate the superconducting gap within the spin- 
fluctuation- mediated pairing mechanism. The gap is basically s±, where the gap changes its 
sign between electron and hole Fermi surfaces, but three-dimensional nodal structures appear 
in the largely warped hole Fermi surface having a strong Z 2 / XZ /Y Z orbital character. The 
present result, together with our previous study of the 1111 systems, explains the strong ma- 
terial dependence of the superconducting gap in iron pnictides. 
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It is now becoming clearer that the superconducting 
gap of iron pnictides 1 has nonuniversal forms. A num- 
ber of experiments have suggested the presence of a fully 
open gap, 2,3 but for LaFePO in particular, experiments 
have suggested the presence of line nodes in the gap. 4-6 
In a previous study, motivated by an experimental ob- 
servation by Lee et al, we have mainly focused on the 
1111 systems, and pointed out that the band structure 
around the wave vector (n, tt) in the unfolded Brillouin 
zone, and thus superconductivity, is sensitive to the pnic- 
togen height measured from the iron plane. 8 Using a five- 
orbital model, 9 we have shown that when pnictogen is at 
high positions, the X 2 — Y 2 band around (tt, tt) gives 
a hole Fermi surface, and in this case, the spin fluctua- 
tions arising from the interaction between electron and 
hole Fermi surfaces give rise to a fully gapped s±-wave 
pairing. 9-12 On the other hand, as the pnictogen is low- 
ered, the Z 2 band around (tt, tt) rises up above the Fermi 
level, while the X 2 ~ Y 2 band sinks below. This results 
in the presence of line nodes in the superconducting gap 
either in the nodal s±-wave form, where the gap nodes 
are on the electron Fermi surfaces, or in the d- wave form 
where the nodes are on the hole Fermi surfaces. 9, 13 Fluc- 
tuation exchange (FLEX) 14 and functional renormaliza- 
tion group 15 ' 16 studies also found a similar tendency. The 
presence of line nodes in the superconducting gap results 
in a lower T c , and this explains the experimental obser- 
vation of line nodes in low-T c LaFePO. 4-6 

However, if we focus on the 122 materials, there are 
now several experiments that do not fit into this view. 
In the isovalent doping system BaFe2(Asi_ 2; P 2 ;)2, 17 ' 18 a 
number of experiments suggested the presence of line 
nodes in the superconducting gap, 19,20 but T c is rela- 
tively high (maximum T c of about 30 K), 18 in contrast 
to LaFePO. Moreover, the X 2 — Y 2 hole Fermi surface 
is found to be present in this material theoretically (as 
seen below), and this also seems to be supported by re- 
sults of an angle resolved photoemission (ARPES) ex- 
periment. 21 Another interesting observation in the 122 
materials is the possible presence of nodes in the super- 



conducting gap of KFe2As2- 22-26 In this material, the 
electron Fermi surface is barely present, 2 so it is likely 
that the gap nodes are on the hole Fermi surface. Here 
also, the X 2 — Y 2 hole Fermi surface should be present 
because there is a large number of holes. These exper- 
imental observations for the 122 materials do not seem 
to fit into the view that nodal pairing occurs when the 
X 2 — Y 2 hole Fermi surface disappears by reducing the 
pnictogen height. As for other origins of the presence of 
line nodes in the gap, the Coulomb avoidance, 28 as well 
as the competition of spin and orbital fluctuations, 29 ' 30 
has been considered. 

In the present study, we consider another possible ori- 
gin of the gap nodes, which is peculiar to the 122 mate- 
rials. We construct a three-dimensional ten-orbital model 
of BaFe2(Asi™ K P ;c )2 by first-principles calculation using 
maximally localized Wannier orbitals, 31 and apply ran- 
dom phase approximation (RPA) to obtain the spin sus- 
ceptibility and superconducting gap function. For the 
122 materials, the Brillouin zone unfolding procedure 9 
that adopts the reduced unit cell (with one iron) can- 
not be performed strictly 32 there fore, by adopting the 
ten-orbital model that uses the original unit cell (with 
two irons) of the body-centered tetragonal lattice struc- 
ture, we fully take into account the peculiar features of 
the 122 band structure not present in the 1111 materi- 
als. The superconducting gap is basically s±, where the 
gap changes its sign between electron and hole Fermi 
surfaces, but when the hole Fermi surface around the Z 
point having a strong 3Z 2 — R 2 (Z 2 ) orbital character be- 
comes large by isovalent doping, the superconducting gap 
on that Fermi surface exhibits three-dimensional nodal 
structures. This kind of node is peculiar to the 122 sys- 
tems, as also found in a five-orbital study of BaFe2As2 
(but with different orbital characters of the Fermi sur- 
face). 33 ' 34 Since the Z 2 orbital does not play an impor- 
tant role in spin-nuctuation-mediated superconductivity, 
this explains why T c in the isovalent doping system is 
relatively high despite the presence of nodes in the su- 
perconducting gap. 
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Fig. 1. (Color online) Band structures of the ten-orbital model of 
(a) BaFe2As2, (b) BaFe2P2, and (c) BaFe2Aso.36Po.64- In (a), 
the original first-principles band calculation is shown by dots. 
The Brillouin zone is shown in the inset. 



In order to construct a realistic model, we first per- 
form first-principles band calculation using the Quan- 
tum Espresso package 35 adopting the experimentally 
determined lattice structure, 36 and then obtain a ten- 
orbital model using maximally localized Wannier or- 
bitals. 31 There are ten bands mainly originating from the 
five 3d orbitals because there are two iron atoms per unit 
cell in the body-centered tetragonal lattice structure. We 
show in the upper panels of Fig. 1 the band structure of 
the ten-orbital models of BaFe2As2 and BaFe2P2- To ob- 
tain a model of BaFe2(Asi_ :c P a ;)2, we first obtain a ten- 
orbital tight-binding model of "hypothetical" BaFe2As2 
and BaFe2P2 having the experimentally determined lat- 
tice structure of BaFe2(Asi_ :E P ;r )2 at each x. We as- 
sume that making a linear combination of the two sets 
of tight-binding parameters (hopping integrals and on- 
site energies) and mixing them with a ratio of 1 — x : x 
will give a good approximation of the band structure of 
BaFe2(Asi_ a; P 2 ;)2. The band structure of this model with 
x = 0.64 is shown in the lower panel of Fig. 1. 

In Fig. 2, we show the Fermi surfaces of 
BaFe2(Aso.36Po.64)2 for each orbital character, where 
the thickness represents the strength of the character. 
The 122 materials share common features with the 1111 
materials in that they have three hole (around T-Z) and 
two electron Fermi (around X-P) surfaces. There are, 
however, some differences. One is that the portion of 
the Fermi surface around the Z point having a strong 
Z 2 character is continuously connected to the XZ/YZ 
portion of the Fermi surface around the T point in 122, 
while in 1111 the Z 2 Fermi surface, when present, is an 
isolated three-dimensional pocket. We will call this hole 
Fermi surface with a mixed Z 2 and XZ/YZ orbital 
character a±. Another difference of 122 from 1111 is 
that the hole Fermi surface having a strong X 2 — Y 2 
orbital character coexists with the Z 2 Fermi surface; in 
1111, either the three-dimensional Z 2 or the cylindrical 
X 2 — Y 2 (7) Fermi surfaces exist depending on the 
pnictogen height. 8 ' 37-39 The hole Fermi surface with an 
X 2 — Y 2 character will be called 7, as in our study of 
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Fig. 2. (Color online) Upper panels: Fermi surfaces of the ten- 
orbital model of BaFc2(Aso.36Po.64)2- Horizontal (vertical) cuts 
are shown in the top (middle) panels. The vertical cuts are pre- 
sented along with the strength of the orbital character. The 
X point corresponds to the wavevector (ir, 0) in the unfolded 
Brillouin zone. Bottom panels: evolution of the Fermi surface 
upon increasing the phosphoruss concentration x. The pnicto- 
gen height hp n is also presented. 



llll. 8 There is another hole Fermi surface having an 
XZ/YZ character, which we will call a^. 

Replacing As by P does not alter the band filling (iso- 
valent doping) , and the main effect is to reduce the pnic- 
togen height measured from the iron planes. In the bot- 
tom panels of Fig. 2, we show the vertical cut of the 
Fermi surface for various P contents. As in 1111, lower- 
ing the height leads to a higher Z 2 orbital level, resulting 
in a larger ai Fermi surface around the Z point. 18 Such 
a strong warping of the Fermi surface has been directly 
observed in a recent ARPES experiment. 40 At the same 
time, the X 2 — Y 2 level is lowered, so that the 7 Fermi 
surface shrinks, but it exists up to a P content of about 
x ~ 0.7 in the present calculation. 

We now move on to the RPA calculation. We mainly 
concentrate on BaFe2(Aso.36Po.64)2 here, while we also 
comment on the calculation results for other x values. As 
for the electron-electron interaction, we consider the in- 
traorbital U, the interorbital V, Hund's coupling J, and 
the pair hopping interaction J'. We consider the orbital- 
dependent interactions obtained in ref. 41 . We apply RPA 
in this model, 8 ' 9 and obtain the spin and charge suscep- 
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Fig. 3. (Color online) Largest eigenvalue of the spin susceptibility 
matrix for the ten-orbital model of BaFe2(Aso.36Po.64)2- The X 
point corresponds to the (it, 0) point in the unfolded Brillouin 
zone. 



tibility matrices. From these, we obtain the pairing inter- 
action, which are plugged into the linearized Eliashberg 
equation. In the RPA calculation (where the self-energy 
correction is neglected) , realistic values of the interaction 
results in very large spin fluctuations, so we multiply all 
the electron-electron interaction by a factor /. 8 In this 
way, we keep the relative strength of the interactions be- 
tween different orbitals to be the same as those obtained 
from the first-principles calculation. In the following, we 
will present the eigenfunction of the Eliashberg equation 
in the band representation, and call them the "(super- 
conducting) gap" . In the actual calculation results shown 
below, we take 16 x 16 x 16 fc-point meshes, 128 Matsub- 
ara frequencies, T = 0.07 eV, and the interaction reduc- 
ing ratio / = 0.55. 

In Fig. 3, we show a plot of the eigenvalue of the spin 
susceptibility matrix. It has a peak at the X point (cor- 
responding to the wave vector (it, 0) in the unfolded Bril- 
louin zone, although the unfolding cannot be strictly per- 
formed), which originates from the interaction between 
electron and hole Fermi surfaces. As shown in our earlier 
study, 8 the main interaction that induces the spin fluc- 
tuation comes from the X 2 — Y 2 and XZ/YZ portions 
of the Fermi surface, namely, the intraorbital interaction 
within these orbitals. Performing a similar calculation for 
other values of x, we find that the spin susceptibility is 
suppressed monotonically upon increasing the phospho- 
russ content, as expected from the fact that the 7 Fermi 
surface shrinks. 8 

We show the superconducting gap in Fig. 4. The elec- 
tron and hole Fermi surfaces have different signs of the 
gap for most portions, namely, the gap is basically s±, 
originating from the spin fluctuation mentioned above. 
The gap on the electron Fermi surfaces is fully open since 
the X 2 — Y 2 -f Fermi surface is present, as is the case for 
the 1111 systems with high pnictogen positions. 8 On the 
other hand, as the volume of the ot\ Fermi surface around 
the Z point grows upon lowering the pnictogen height 
by isovalent doping, a three-dimensional sign change of 
the gap takes place within this Fermi surface, as shown 
in Fig. 4. The gap function on the a± Fermi surface is 
schematically shown at the bottom of Fig. 4. This sign 
change can be considered to be due to the repulsive in- 
traband interaction within the a± Fermi surface (shown 
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Fig. 4. (Color online) Contour plots of the gap function on the 
five Fermi surfaces for five horizontal cuts around T— Z (hole) or 
X— P— X (electron). The solid lines represent the Fermi surfaces, 
while the dashed lines are the nodes of the gap. The Brillouin 
zone is shown along with the regions where the contour plots of 
the gap are presented. In the bottom, a schematic view of the gap 
function on the ai=Z 2 / XZ /Y Z hole Fermi surface is shown. 
The arrow bridging the Fermi surface indicates the intraband 
interaction that induces the sign change of the gap. 



by the arrow in the schematic figure), which becomes 
more effective as the volume of the Fermi surface around 
the Z point grows. In fact, we find that the plus-sign re- 
gion of the gap on the a\ Fermi surface tends to shrink 
for smaller values of x. A similar three-dimensional sign 
change of the superconducting gap has been found in 
a five orbital study of Bal22 in ref. 33 , although in the 
model of ref. 33 . the a\ Fermi surface around the Z point 
has an X 2 — Y 2 character (or xy in the notation of ref. 33 
.) rather than Z 2 . 
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Since the three-dimensional nodes in the present study 
occur at portions of the Fermi surface having a strong Z 2 
character, the presence of the nodes does not strongly re- 
duce T c . This is in contrast to that in the case of LaFePO, 
where nodes of the gap, according to a similar RPA and 
other theoretical studies, enters in the X 2 — Y 2 and/or 
XZ/YZ Fermi surfaces, 8 ' 14-16,42 which play an impor- 
tant role in spin-fluctuation-mcdiatcd superconductivity. 
The present view explains why T c is relatively high in 
BaFe 2 (Asi_ :E P a ;)2 despite the presence of nodes in the 
superconducting gap. 

To summarize, we have obtained a three-dimensional 
ten-orbital model for the isovalent doping material 
BaFe2(Asi_ :E P a ;)2, and applied RPA to obtain the su- 
perconducting gap originating from spin-fluctuation- 
mediated pairing. The obtained gap is basically s±, and 
we find no gap nodes on the electron Fermi surface be- 
cause the X 2 — Y 2 7 Fermi surface persists up to a high P 
content. On the other hand, there are three-dimensional 
nodes on the Z 2 /X Z/YZ portion of the hole Fermi sur- 
face. We presume that these three-dimensional nodes in 
the Z 2 /XZ/YZ hole Fermi surface are responsible for 
the experimentally observed nodal behavior of the super- 
conducting gap of BaFe2(Asi_ a; P a ;)2, since the presence 
of nodes in this material does not seem to strongly af- 
fect T c . The relation between the present view and other 
experimental results will be of great interest. A compari- 
son between results of neutron scattering experiment and 
theoretical analysis based on the present gap structure 
is now in progress. 43 It is also necessary to understand 
how such a gap structure on the hole Fermi surface can 
be consistent with results of the ARPES experiment, 21 
or the apparently "nodeless" behavior observed in a spe- 
cific heat experiment. 44 We note here that our study 
does not completely rule out the possibility of nodes on 
the electron Fermi surface, since correlation effects that 
are not taken into account in the present "band calcula- 
tion+RPA" study may give rise to gap nodes on the elec- 
tron Fermi surface (even when the X 2 — Y 2 hole Fermi 
surface is present). 14 Also, there is a possibility that in 
the efeciron-doped 122 materials, where the 7 Fermi sur- 
face is less effective, the nodes may be on the electron 
Fermi surface, as suggested by a recent study. 32 In any 
case, the tendency of three-dimensional nodes entering 
the Z 2 /XZ/YZ hole Fermi surface of the 122 materi- 
als should remain as long as this Fermi surface is suffi- 
ciently large and the repulsive intraband interaction is 
effective. In this context, a "heavily hole-doped" mate- 
rial, KFe2As2, mentioned as another nodal superconduc- 
tor in the introductory part, is also of great interest. In 
fact, in our preliminary study of this material, we find 
a similar tendency of three-dimensional gap nodes enter- 
ing the Z 2 /XZ/YZ hole Fermi surface. In this study, we 
use a model that quantitatively reproduces the spin fluc- 
tuation modes observed in a recent neutron scattering 
experiment. 45 The details will be published elsewhere. 
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